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Abstract 

*ys The use of resistive technologies to MPGD detectors is taking advantage for many new applications, including high rate and 

t— i energetic particle flux scenarios. The recent use of these technologies in large area detectors makes necessary to understand and 

<^i characterize the response of this type of detectors in order to optimize or constrain the parameters used in its production, material 

-J resistivity, strip width, or layer thickness. The values to be chosen will depend on the environmental conditions in which the detector 

^h will be placed, and the requirements in time resolution and gain, improving the detector performance for each given application. We 
^TJ" 1 present two different methods to calculate the propagation of charge diffusion through different resistive topologies; one is based 

"^h on a FEM of solving the telegraph equation in our particular strip detector scheme, the other is based on a semi-analytical approach 

(**-« of charge diffusion and is used to determine the charge evolution in a resistive layer. 



i Keywords: micromegas, charge diffusion 
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1. Introduction 

One of the most prominent problems associated with gas- 
filled proportional chambers is the sparking induced by heavily 
ionizing particles producing large deposits. Amplified by the 
avalanche process they could reach a critical charge density, 
related to the Raether's limit [1|, and could evolve into a dis- 
charge. As a result, a large fraction of the stored charge defin- 
ing the amplification field is lost. The presence of discharges 
limits the high rate operation of the detector due to the required 
power supply recovery time, reduces the mean life of the detec- 
tor and risks the damage of the readout electronics due to the 
high currents reached by these spark processes. 

To avoid any damage of the electronics, most of the detectors 
are using an additional protecting circuit, which interfaces the 
readout strips or pads with the front-end electronics. Recently 
new developments are undertaken to further improve the spark 
protection to operate at higher rates. 

A promising approach is employing a resistive material on 
top of the anode plane [2|. When a streamer develops the in- 
creasing number of charges is deposited at the high resistivity 
material. The anode potential generated by these charges re- 
duces the amplification field, quenching the streamer resulting 
from field loss. A required condition is that the anode potential 
must last for a period of time long enough such that the charges 
present in the streamer have been completely evacuated, which 
introduces the need for high resistive values. The resistive in- 
terface limits the charge transfer, protecting the electronics, and 
avoids the total discharge of the mesh reducing the dead time of 
the detector. 

This new technology needs to be characterized and simulated 
to exploit its possibilities, using the appropriate resistive and 



capacitive values depending on the needs of each particular ap- 
plication. 

Micromegas detectors have evolved recently using this kind 
of protection. Two different topologies are typically imple- 
mented, a fully covering resistive layer and micromegas detec- 
tors with resistive strips. We present the description of these 
new approaches focusing on recent developments. Then, we 
will introduce a method to calculate the charge dispersion on 
the resistive material in each topology. 

1.1. Resistive strip topology 

For large Micromegas detectors, a novel protection scheme 
has been developed by a group working on hard radiation detec- 
tors at CERN 1 3 , 4 1 . According to this scheme, we are using re- 
sistive strips above the readout electrode to overcome the spark 
problem. This technique consists in covering the micromegas 
readout |5 , 6| with a thin insulating layer where resistive strips 
are plotted on one favored axis (see Figure fTJ. 
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Figure 1 : Schema of the resistive micromegas detectors 
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By adding a drift electrode, a MM detector operates as usu- 
ally in the proportional avalanche mode inducing signals on the 
resistive strips plane. The resistive strip is grounded through an 
embedded high impedance resistor. The typical strip resistivity 
is in the order of 100 MQ/cm, and 200 MQ. for the boundary 
resistor which drives the strip to ground. 

1.2. Resistive layer topology 

An approach of resistive layer topology is called Piggyback 
resistive Micromegas Q where a thin resistive layer is de- 
posited on an insulator. This element is used as substrate for the 
Micromegas which could be fabricated by using the bulk tech- 
nology [8|. A schematic view of the new structure is shown 
in figure [2] where we can see the various elements: the Mi- 
cromegas structure with the pillars; the resistive layer; the in- 
sulating substrate and the readout board. The thickness of the 
resistive layer is of the order of 10 /mi. 
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Figure 2: Schema of the Piggyback detectors 



A particularity of this new detector technique is the fact that 
the readout electronics are independent (and thus exchangeable) 
from the micromegas structure. The elements of the structure 
are optimized in such a way that the electronic signal is not 
lost through the resistive layer but is propagated to a separate 
plane, carrying read-out pads or strips, by capacitive coupling. 
In order to optimize the induced signal by capacitive coupling 
the thickness of the insulating ceramic (b) must be kept small 
to satisfy the relation 



h <c h - 



the amplification gap (ti) is about 100 fim and it is achieved 
between a woven stainless steel mesh connected to the high 
voltage and a resistive anode. e\ is the dielectric constant of 
medium 1 (gas) and e-i is the dielectric constant of medium 2 
(insulator). Because the insulator plays at the same time the 
role of vessel of the detector the thickness should be kept rea- 
sonable (several hundred //m). In order to satisfy ([T]i the mate- 
rial should have a dielectric constant as high as possible. Good 
candidates are ceramic insulators having large dielectric con- 
stants (»10). The value for the amplification gap is the stan- 
dard 128yum, while the ceramic thickness should be thicker than 
300/mi to assure the leak tightness of the detector. The resis- 
tive layer is made of ruthenium oxide (RUO2). This material is 
extensively used for coating ceramics at high temperature for 
the preparation of resistors or integrated circuits, and provides 
good uniformity and a wide range of resistivities, between 1- 
lOOMn/D. 



2. Simulation of charge diffusion through a resistive strip 

A simplified approach of a single resistive strip is used in 
order to calculate the induced signal at the micromegas usual 
read-out and the charge diffusion through the resistive strip. 
The model here described is based on the equivalent electric 
circuit introduced by [9]. This equivalent circuit describes in 
general a standard micromegas detector covered by a thin insu- 
lating layer where resistive strips are placed facing the standard 
metallic strips, creating an inter-capacitance and therefore an 
electric transmission line. 

In order to calculate the charge diffusion, the resistive strip is 
subdivided in a set of differential elements which final element 
is grounded through a boundary resistor Rb- The resistive strip 
has a linear resistivity R A and is electrically coupled to the stan- 
dard read-out, defined by V c , through the mentioned insulat- 
ing layer which introduces the strip linear capacitance C\. The 
standard read-out is connected in this model to a low impedance 
value R s t r i P , and the residual capacitance of the PCB board C pc h 
(see figure[3]l. 
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Figure 3: Simplified RC-model of a resistive strip. 



The potential on the strip is driven by the following equation 



d 2 V(x,t) 
dx 2 



= C A R, 



d(V(x,t)-v c (t)) 
at 



(2) 



The evolution of the potential and/or charge in the resistive 
strip is bounded by the external circuit elements attached to the 
(1) strip electrode, R slr i P and C pc b, and the fixed resistive value, 
Rb, attached to the resistive strip. The first elements {ft strip an d 



C pc b) allow determination of the strip voltage V c as a function of 
the resistive strip potential distribution, by integrating the total 
equivalent current induced at the capacitive surface ic A along 
the resistive strip length x L , 



i Cl =C A £ L §- t {V{x,t)-V c {t))dx = 

dV c (t) , V c (t) 



-P cb dt 



(3) 



tfw 



The resistive value Rb provides the boundary condition on 
the voltage V„{t) at the beginning of the strip, and its value is 
related with the voltage gradient at the resistive strip end, 



V {t) 



R b 8V{x,t)\ 



R, 



dx 



(4) 



The relation pi describes the potential propagation at the re- 
sistive strip. However, in order to observe such propagation we 



need a perturbation which is introduced as a time dependent 
current source spatially distributed and described by the change 
on linear charge density p(x, f) along the resistive strip. This 
current is introduced inside relation (ph by taking into account 
the expected voltage increase at the resistive nodes, 



d 2 V(x,t) 
8x 2 



CaRa 



d(V(x,t) 



dt 



V c (0) „ dp(x,i) 

+R *~dT 



(5) 



the introduced charge density pretends to emulate the effect of 
charge induced by a physical avalanche, which would be the 
connection between our detector read-out and the physical sig- 
nal. 

The solution of the coupled relations ([3]) and p) for a given 
input charge density p(x, f) will provide a full description of the 
signal propagation within the frame of this model (a detailed 
solution can be found at 1 10 1). 

This model allows to solve particular scenarios and obtain the 
transition state of the strip receiving a charge stimulus. Figure HI 
shows the readout signal at the resistive boundary resistor and 
its shape dependency with a current applied at different strip 
positions for a total strip length of 10 cm. 
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Figure 4: Signal observed at the boundary resistor for currents applied at dif- 
ferent resistive strip positions. 



This method can be also used to obtain the potential profile 
in the resistive strip for different rate configurations. Figure B] 
shows the final state potential distribution resulting from an ho- 
mogeneous rate distribution. 
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Figure 5: Final potential distribution along the resistive strip for an homoge- 
neous current distribution. 



3. Simulation of charge diffusion effects over a resistive 
plate 

In this section we show the effect on the gain due to the ac- 
cumulation of charge over time due to slow charge diffusion in 
the illuminated region for a resistive layer topology. The charge 
diffusion relation for a given charge Q, which follows a Gaus- 
sian distribution of width w, was provided in ifTTI 



p(r, t) 



Q 



2n{2ht + w 2 ) 



exp 



2(2fo + w 2 ) 



h= l/RC 



where the parameter h it is related to the surface resistivity R 
and surface capacitance C. The accumulated charge will be 
translated into a surface potential given by the capacitance of 
the ceramic and thus, in a reduction of the amplification field. 
In order to study such reduction a narrow region around the 
Gaussian distributed charge has been chosen. The evolution on 
time of the total charge inside this region, defined by the radius 
R , it is integrated from the previous expression and is described 
by the following relation 



p(Ro,t) = Q 



1 - exp 



Ri 



2(2ht + w 2 ) 



We can extrapolate this relation to the case of a continuous 
current flow by assuming this relation is valid for a differential 
element of charge. Integrating these differential contributions 
over time will allow us to calculate the accumulated charge at 
any time, 



Jo 



Q(R ,to) = fo° 



dpQU) 



1 - exp 



( 



dt 



K 



2{2h{t a -t)+w 2 ) 



dt = 



(6) 



where we should take into account that these differential charge 
contributions will be related to the local amplification field, 
modified by the anode potential V a at the resistive plane at each 
time interval. 

The integral will be calculated by adding the charge contri- 
butions at each time interval n, and describing the gain in terms 
of the relative mesh versus anode potential, 



Sq„ = g(V a (n6t))N e q e r 



VJnSt) = qJnRlC 



where N e is the number of electrons produced in each interac- 
tion, r is the interaction rate, g is the gain as a function of the 
anode potential V a , The gain as a function of the mesh poten- 
tial is obtained experimentally by illuminating the Piggyback at 
low rates (see figure 16k expression that will be used afterwards 
with the corrected potential given by V a , 

The charge is integrated considering the effect the accumu- 
lated charge will have in the gain at each time step. The sim- 
ulation of each time step requires to know previously all the 
charges at each past time interval in order to apply the diffu- 
sion relation to them, this makes the computation of each time 
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Figure 6: Absolute gain measured (blue points) versus applied mesh voltage in 
Ar+10% C02, and the fitting result (red-filled curve). 



interval more expensive and the time required has a non-linear 
growth. The simulation of the first half an hour with reasonable 
accuracy can still be performed within a few hours. 

A set of simulations was launched at different scanning rates 
for different RC values in order to observe the effect on the 
gain drop. Figure|7]shows the final gain value obtained after an 
exposure of 5400 seconds, several simulations were performed 
at different rates in order to show the result as a function of the 
interaction flux. Different resistivity and capacitive values are 
tested from the range of 10MQ/D to lOGfi/D, at 1 pF/mm 2 . 
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Figure 8: Experimental gain drop dependency with the applied rate for different 
applied mesh voltages, V,„ . 



We showed two methods which allows us to calculate the 
charge diffusion on the resistive material and predict the effects 
on the gain in each topology. The simulated effects are exper- 
imentally observed for the case of Piggyback detectors. The 
knowledge of the relation between the resistive values and the 
expected gain drop as a function of the rate will allow us to de- 
termine the optimum parameters to be used depending on the 
final application requirements. 
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Figure 7: Relative gain versus the particle flux for different resistivity values 
at lpF/mm 2 . 

As lowest is the resistivity of the material, the longest is the 
fiat gain region, and higher rates could be achieved without ex- 
pecting an effect on the gain. 

Experimental measurements performed with an X-ray gener- 
ator show a behavior similar to the simulated one (see FigurelH). 
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4. Summary 

We report on two recent representative developments on re- 
sistive micromegas technology based on resistive layer and re- 
sistive strips topologies. These techniques profit from the resis- 
tive coating by increasing the mean life of the detector, reduc- 
ing the sparks damage of the electronics and improving the rate 
capability by reducing the full discharge of the cathode. 



